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Abstract: Hibiscus aridicola (Malvaceae) is an endangered ornamental shrub endemic to 
the dry-hot valleys of Jinsha River in southwest China. Only four natural populations of 
H. aridicola exist in the wild according to our field investigation. It can be inferred that 
H. aridicola is facing a very high risk of extinction in the wild and an urgent conservation 
strategy is required. By using a modified biotin-streptavidin capture method, a total of 
40 microsatellite markers were developed and characterized in H. aridicola for the first 
time. Polymorphisms were evaluated in 39 individuals from four natural populations. 
Fifteen of the markers showed polymorphisms with two to six alleles per locus; the 
observed heterozygosity ranged from 0.19 to 0.72. These microsatellite loci would be 
useful tools for population genetics studies on H. aridicola and other con-generic species 
which are important to the conservation and development of endangered species. 

Keywords: Hibiscus aridicola; endangered plant; SSR markers; population structure; 
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1. Introduction 

Genus Hibiscus (Malvaceae) includes at least 250 species distributed in tropical and subtropical 
regions [1]. H. aridicola is endemic to the dry and hot valleys of the Jinsha River (upper reaches of the 
Yangtse River) in Southwest China. It is an ornamental deciduous shrub with large flowers that vary 
from white to whitish purple. Although the species was described in 1927 [2], and literature or relative 
research references are limited. H. aridicola was widely distributed along the valleys of Jinsha River 
from 1300 m to 2100 m during the 1950s [2-4]. However, it currently has been evaluated as 
Endangered [EN B2ab (if)] (following IUCN Red List Categories and Criteria) with less than five 
locations and a continuing decline of area of occupancy [5]. Our recent field surveys have further 
confirmed the threatened status of the species and only four natural populations ofH. aridicola exist in 
the wild. Our surveys have also discovered that habitats of H aridicola have been greatly degraded 
and the populations are fragmented and isolated from each other. Population genetic studies have 
predicted that fragmentation will lead to a loss of genetic diversity due to inbreeding, population 
isolation and restricted gene flow and small effective population sizes and that may lead to a decline in 
fitness or even, ultimately, extinction [6,7]. As an endangered and endemic ornamental shrub, 
H. aridicola is facing a very high risk of extinction in the wild, and its effective and long-term 
conservation is urgently needed. In our study, 40 micro satellite markers were developed and 
characterized to investigate the genetic diversity and population structure using the fast isolation by 
Amplified Fragment Length Polymorphism (AFLP) of sequences containing repeats (FIASCO) [8]. 

2. Results and Discussion 

A total of 441 positive clones were captured, among these 232 clones (53%) were found to contain 
simple sequence repeat (SSR). Finally, 84 sequences contained SSR loci were selected for primer 
design. 40 micro satellite loci successfully amplified in H. aridicola. We tested the degree of 
polymorphism in four populations of H. aridicola for 40 microsatellite loci and 15 out of them are 
polymorphic amplification, other 25 microsatellite loci were monomorphic by the result of 
polyacrylamide gel (Table 1). 



Table 1. Primer sequences and characteristics of 40 microsatellite loci successfully 
amplified in 39 Hibiscus aridicola individuals. 



Locus 


Primer Sequence (5'-3') 


Repeat 
motif 


NA 


Size 
(bp) 


(°C) 


Gene Bank 
Accession No. 


HA-1 


F: TTGAACATAAACAAGCGG 
R: AAAACAAGTTGGGGAGG 


(GT) 14 


2 


180-200 


56 


JN167557 


HA-2 


F: CTGAATGCCAGAATGACT 
R: CAGGCGAAAGAGGAAGAT 


(AG) 21 


5 


410-445 


60 


JN167558 


HA-3 


F: ATCATTATCATCTTCGTTTC 
R: AAGGGACCAAAGTCTCAA 


(CT) 15 


2 


141-154 


59 


JN167559 


HA-4 


F: CACCAAATCCTGGAGAAG 
R: GCAAACGAGAATAATCAAAA 


(AC)„ 


2 


137-150 


60 


JN167560 
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Table 1. Cont. 



Locus 


Primer Sequence (5'-3') 


Repeat 
motif 


N A 


Size 
(bp) 


(°C) 


Gene Bank 
Accession No. 


HA-5 


F: GCGTGGATGTTCTTTCTT 
R: TCGAACCCTATGGATGTA 


(GA) 8 


3 


155-175 


61 


JN167561 


HA -6 


F: GAACAAGCCTGTCACTAA 
R: CACAAACCGATTTACGAT 


(CT) 18 


2 


130-178 


60 


JN167562 


HA-7 


F: CAGCAGTTAGAGCAGGAGGT 
R: TTCGGACATGAGTATGGGAT 


(TG) 8 


3 


208-215 


59 


JN167563 


HA-8 


F: CACTTCCACGAAGCTCTTAC 
R: GGAGATAAACAGAAAAGGGTA 


(CTTCT) 3 


2 


230-242 


59 


JN167564 


HA-9 


F: TATGGGTTTAGTGCCTGTAT 
R: TAGGTTGCTTGAATCTTTTC 


(AC), 


2 


322-341 


59 


JN167565 


HA- 10 


F: CCCAAACCTCTATCATCT 
R: ATATCCCTTAGTTCTGCT 


(GT)„ 


2 


197-204 


59 


JN167566 


HA-11 


F: CACCAAATCCTGGAGAAGTA 
R: GGCAAACGAGAATAATCAAA 


(AC), 


3 


98-109 


58 


JN167567 


HA- 12 


F: AAGGAGAAGCCAAGGTGAAA 
R: GACAAACCCACATACAGGAA 


(GAA) 5 


2 


119-132 


60 


JN167568 


HA- 13 


F: ACTTTTATCGTATAGACCAG 
R: G AAC AC CTTTATTTC AGTGT 


(CTT) 15 


2 


110-118 


59 


JN167569 


HA- 14 


F: GAAATGGCAAGGTTTTAGAT 
R: CTCAACTTTTGTGATGTGGC 


(GAA)i 0 


3 


144-153 


59 


JN167570 


HA- 15 


F: CAGCCACAATCCTCCTAACT 
R: GAAGGGTAACTTGTTTCACG 


(TG)„ 


3 


313-345 


60 


JN167571 


HA- 16 


F: TTGAGATTTG AC CTGGAA 
R: AC ATTGGC GAAGATAC AC 


(CT) 20 


. 


237 


57 


JN167572 


HA- 17 


F: TATTTCCCTGTCCCTGTT 
R: GACCTTTTCGTCTTTTGG 


(CT) 14 




74 


54 


JN167573 


HA- 18 


F: CACCCAAGCATGATAAAA 
R: AGAATGAAAGAAAATGGC 


(AC) 14 




164 


58 


JN167574 


HA- 19 


F: ACCACCAGAAAGCAAACA 
R: GATGACTAATGGGAAAGAA 


(CT)i2 


i 


122 


58 


JN167575 


HA-20- 


F: TCGTGATGGGAACAGATA 
R: TGAAATACTCATGGGAATG 


(GA) 13 


i 


141 


57 


JN167576 


HA-21 


F: AGAAAATCCCAATCTCAA 
R: CTAGCCAGAAACAACGAG 


(CT) 16 


. 


201 


61 


JN167577 


HA-22 


F: ACTGGTAACATCCCTGAC 
R: GAAACTGCTGGAAATCTA 


(AC), 


i 


107 


60 


JN167578 


HA-23 


F: AGCATCCGATCCTTATCT 
R: TATCAGCGACTCCTCCAC 


(CT) 16 




158 


59 


JN167579 


HA-24 


F: AGTCATCGGAGAAATAGAG 
R: ATAACCAAGGAGGAAACA 


(CAT) 9 




442 


58 


JN167580 


HA-25 


F: AAACTGCGAAATCCTCAT 
R: AGTAAACACTGCCTCCAT 


(CT) 16 (CA) 8 




100 


60 


JN167581 
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Table 1. Cont. 



Locus 


Primer Sequence (5'-3') 


Repeat 
motif 


N A 


Size 
(bp) 


(°C) 


Gene Bank 
Accession No. 


HA-26 


F: CCTCCGTGGTAACTCCTT 
R: TGATGAAATATGGCTTGG 


(CA)„ 


1 


137 


59 


JN167582 


HA-27 


F: TGAATTTCTTTTCTTC CTTTAC 
R: CAACTATCATCTTGTCGTGC 


(TG) 10 


1 


207 


58 


JN167583 


HA-28 


F: ATAGATG A AC C AGG AA AT 
R: CTGAAGATAAAGAAAGCA 


(AC) 10 


1 


102 


53 


JN167584 


HA-29 


F: ATACGACAGATGCGGAAGTG 
R: TTAGTTACGGGAAC CGAAGG 


(GAA) 3 


1 


196 


60 


JN167585 


HA-30 


F: TTGCTCACTTGAAAACATTA 
R: GAAAACGACACGATCACTCT 


(AC) 10 


1 


72 


58 


JN167586 


HA-31 


F: GGAAAGTGGCTGACTGGTAG 
R: C GAC ATCGGTG AGGTTGGTT 


(AC), 


1 


296 


60 


JN167587 


HA-32 


F: ACGGAAATGCTCAAACCCTC 
R: AAATGATTACCGCCGACAAC 


(CTT) 3 


1 


121 


58 


JN167588 


HA-33 


F: GCTC AGGTAAACC C ATAA 
R: GCTCGTCGTACATACACTT 


(CT) 7 


1 


300 


55 


JN167589 


HA-34 


F: TACTGTCCAATGAATGCCTT 

R: AACCTGAACTATAAAGTAAACTGC 


(GT) 14 


1 


155 


60 


JN167590 


HA-35 


F: GTATGTTGCTATCCCCTAT 
R: CAAACCAAACAACACTAA 


(TG) 8 


1 


139 


60 


JN167591 


HA-36 


F: GAAAGGAATTGTACGTGGCA 
R: TATGGCTTGGGATTGGTTTT 


(CA) 8 


1 


264 


60 


JN167592 


HA-37 


F: TAAGATGGTATTGGAAGGG 
R: AGGGAGCATAAAAGTGGT 


(GT)„ 


1 


345 


60 


JN167593 


HA-38 


F: ATCATCGGCAGCGACTAG 
R: CAGCAAGGACATCAGGGT 


(GA)i 3 




203 


59 


JN167594 


HA-39 


F: TCCAAGATACTGCCATAC 
R: GGTTCTACAGGTACATGC 


(AC) 14 




108 


60 


JN167595 


HA-40 


F: GAAGAGC GAC AGAAAATG 
R: GGAAAATAAACAAGGGTAAA 


(AC) 14 




103 


59 


JN167596 



T. d , PCR annealing temperature; N A , number of alleles; Size range, size of alleles. 



The number of alleles ranged from two to six in 39 individuals of the species sampled from the four 
extant natural populations. Value for Ho and He ranged from 0.19 to 0.72 (mean Ho = 0.52) and from 
0.20 to 1.00 (mean// E = 0.62), respectively (Table 2). 

Table 2. Results of 15 polymorphic primers screening in four populations ofH aridicola. 



Population 1 (N = 10) Population 2 (N = 10) Population 3 (N = 10) Population 4 (N = 9) 



Locus 


N A 


Ho 


H E 


N A 


Ho 


H E 


N A 


Ho 


H E 


N A 


Ho 


He 


HA-1 


2 


0.51 


0.70 


2 


0.51 


0.80 


2 


0.48 


0.70 


2 


0.47 


0.67 


HA-2 


5 


0.72 


1.00 


4 


0.49 


0.60 


3 


0.59 


0.60 


4 


0.66 


0.78 


HA-3 


2 


0.19 


0.20 


2 


0.48 


0.10 


2 


0.47 


0.40 


2 


0.52 


0.00 


HA-4 


2 


0.51 


0.80 


2 


0.39 


0.50 


2 


0.48 


0.70 


2 


0.53 


1.00 
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Table 2. Cont. 



Population 1 (N = 10) Population 2 (N = 10) Population 3 (jV = 10) Population 4 (jV = 9) 



Locus 


A J 


T T 

H 0 


TT 

H E 


A T 


TT 

H 0 


TT 

He 


A T 


TT 

H 0 


TT 

H E 


A T 


T J 

Ho 


TT 

He 


HAo 


3 


U.62 


U. /O 


3 


0.5 / 


n a n 
0.4U 


3 


n ga 
U.54 


U. /O 


J 


n a £. 
U.46 


U.D6 


HA-6 


2 


0.48 


0.70 


2 


0.39 


0.50 


2 


n ao 

0.48 


0.70 


2 


0.53 


1.00 


I I A -7 

HA- / 


3 


A A "7 

0.4 / 


U.5U 


L 


C\ A A 

0.44 


U.6U 


3 


0.56 


n fid 


L 


U.J / 


n a a 
0.44 


HA-8 


2 


0.52 


0.70 


3 


0.57 


0.60 


3 


0.65 


0.70 


2 


0.53 


0.78 


HA-9 


2 


0.44 


0.60 


2 


0.52 


0.70 


2 


0.44 


0.60 


2 


0.50 


0.78 


HA- 10 


2 


0.39 


0.50 


2 


0.44 


0.60 


2 


0.39 


0.50 


2 


0.37 


0.44 


HA- 11 


3 


0.51 


0.70 


2 


0.51 


0.60 


3 


0.62 


0.70 


3 


0.69 


0.78 


HA- 12 


2 


0.44 


0.60 


3 


0.66 


1.00 


3 


0.65 


0.90 


3 


0.60 


0.89 


HA- 13 


2 


0.39 


0.30 


3 


0.57 


0.60 


3 


0.42 


0.30 


3 


0.65 


0.56 


HA- 14 


3 


0.42 


0.40 


4 


0.57 


0.30 


3 


0.43 


0.20 


3 


0.57 


0.44 


HA- 15 


3 


0.53 


0.70 


3 


0.61 


0.70 


3 


0.68 


1.00 


3 


0.68 


0.89 



N = population samples size; N A , number of alleles; H 0 , observed heterozygosity; H E , expected 
heterozygosity; Population 1 (27°15'N, 102°53'E); Population 2 (27°44'N, 100°24'E); Population 3 (27°45'N, 
100°18'E); Population 4 (27°29'N, 100°09'E). 



These micro satellite markers developed in our study will be a useful tool for further studies of 
conservation genetics, and will help us understand the genetic structure of H. aridicola, so as to make 
effective conservation strategy for this endangered plant. 

3. Experimental Section 

Total genomic DNA of H. aridicola was extracted from dry leaf tissue of three different individuals 
using a modified 4 x CTAB method [9]. A microsatellite enriched library was conducted by using a 
modified biotin-streptavidin capture method. Firstly, the genomic DNA (about 800-1000 ng) was 
completely digested with Msel restriction enzyme (NEB) into the sequences between 200 bp-800 bp, 
then the digested fragments were ligated to Msel AFLP adaptor following by amplification with 
adaptor-specific primers (5'-GAT GAG TCC TGA GTA). After a PCR the products were hybridized 
with biotinylated probes (AG)is, (AC)i5 and (AAG)io respectively [8]. The PCR products were ligated 
into PGEM-T vector (Promega) after purified, and transferred into E. coli strain DH5a (Tiangen). The 
proportion of products and PGEM-T vector in the total volume should be controlled between 3:1 and 
5:1. Positive clones were picked out and tested using (AG)io/(AC)io/(AAG)7 primers and vector 
primers SP6/T7 respectively to select appropriate fragments which contained SSR. The positive clones 
were captured for sequencing with an ABI PRISM 3730XL SEQUENCER. Sequences contained 
simple sequence repeat were selected for primer design using Primer Premier 5.0 program [10]. 

The designed Primer pairs were tested in 39 wild H. aridicola samples from four populations in 
Southeast China. Herbarium voucher deposited in Kunming Institute of Botany, Chinese Academy of 
Sciences (code ZL00 17-00 19, ZL0108- ZL 0120). The PCR amplification was carried out in a total 
volume of 10 jj,L reaction containing 5 \iL 2 x Taq PCR MasterMix (Tiangen; 0.1 U Taq 
Polymerase/nL, 0.5 mM dNTP each, 20 mM Tris-HCl (Ph8.3), 100 mM KC1, 3 mM MgCl 2 ), 0.5 |xL of 
each primer and 0.5 [iL Genomic DNA. Amplification was carried out in thermocycler (veriti 96 well 
thermal cycler) with a cycling profile 95 °C for 4 min then 35 cycles of 45 s at 94 °C, 45 s at the 
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specific annealing temperature, 45 s at 72 °C and a final extension step of 8 min at 72 °C. The 
amplification products were separated on 8% denaturing polyacrylamide gels and visualized by silver 
staining with a 100 bp extended DNA ladder (Fermentas) as a size standard. 

The data was analyzed by GENEPOP 4.0 [11], which included test of observed heterozygosity (Ho) 
and expected heterozygosity (H E ) for the 15 polymorphic micro satellite loci. 

4. Conclusions 

In summary, 40 micro satellite markers have been specifically developed for H aricicola in this 
study. The high discriminatory power of 15 polymorphic loci suggests that they should be suitable for 
the fine-scale analysis and survey of population structure in scarce populations of the endangered 
H aridicola. These developed and characterized SSR markers for H aridicola would also be useful 
for exploring genetic diversity and genetic structure of other species in Hibiscus. 
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